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We report a precise determination of the B^-B'^ mixing parameter Am^ based on the time 
evolution of same-sign and opposite-sign dilepton yields in T(45') decays. Data were col- 
lected with the Belle detector at KEKB. Using data samples of 29.4 fb~^ recorded at the 
T(45') resonance and 3.0 fb~^ recorded at an energy 60 MeV below the resonance, we measure 
Amd = (0.503 ± 0.008(stat) ± O.OlO(sys)) ps~^. From the same analysis, we also measure the ratio 
of charged and neutral B meson production at the T(45'), f+/fo = 1-01 ± 0.03(stat) ± 0.09(sys), 
and CPT violation parameters in B°-B° mixing, K(cos6') = 0.00 ± 0.12(stat) ± O.Ol(sys) and 
Q{cose) = 0.03 ± O.Ol(stat) ± 0.03(sys). 

PACS numbers: 13.66.Bc, 13.25.Gv, 14.40.Gx 



INTRODUCTION 

The mass difference of the B^-B^ mass eigen states, 
Amd is a fundamental parameter in the B meson system. 
The techniques that have been employed to measure it in 
T(45') B^B^ decays fall into two categories, namely 
inclusive and exclusive. Among the inclusive methods, 
the analysis of events where both B mesons decay into 
a final state that includes a high momentum lepton pro- 
vides the largest event sample and is well suited for 
further high precision measurements of the time evolu- 
tion in the B^B^ system. This system exhibits sensitiv- 
ity not only to Am^, but also to other potentially inter- 
esting phenomena such as CP violation in mixing, the 
decay width difference between the two mass eigenstates 
and possible CPT violation in mixing 

Without the assumption of CPT invariance, the 
flavour and mass eigenstates of the neutral B mesons 
are related by 

\Bh) = p\B°)+q\B') 
\Bl) = p'\B")-q'\B''). (1) 
The coefficients p, q, p' and q' can be expressed in terms 



of the complex parameters 9 and (phy ^ — tan(|)e*'^ and 
±. = cot(§)e*'^. CP is violated if ^ 0, and CPT is 
violated if ^ ^ f The time-dependent decay rates 
are given by [2| 



T(4S) 



,-|At|/r^o. 



I sm0e^"''f [cosh (^Ai) - cos {ArridAt)] (2) 
for same-sign dilepton (SS) events, and 



r 



At|/TgO 



4rB0 
' Ar 



B" X 



T{4:S)^e+i- (At) 

{(1 + |cos6'p)cosh(^At) + (1 - |cos6lp)cos(AmrfAt) 
+ 2SR(cos e) sinh At) - 23(cos 6) sm{AmdAt)} (3) 

for opposite-sign (OS) events. Here, we assume CP 
is conserved in the flavour specific semileptonic decay 
amplitudes of the neutral B mesons and set Ai = 
{X-e+Vi\B°) and Ai = {X+£-Pe\B°) to be equal. Am^ 
and Ar are the differences in the mass and decay width 
between the two mass eigenstates of the neutral B meson, 
r — l/rgo is the average decay width of the two mass 
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eigenstates, At is the proper time difference between the 
two B meson decays and is defined as At = t(£^) — 1{£~) 
for the OS events, while the absolute value is taken for 
SS events. 

In equation 121 CP violation appears as a difference in 
the £^£^ and £~i~ rates, in case of a non-zero It 
does not depend on 9, AT or the fraction of mixed events. 
The last two terms in equation |31 are clearly asymmetric 
in At. The last term will dominate over the second-to- 
last term since Arud 3> AF. In this analysis, we assume 
Ar and CP violating effects are negligibly small 0, Q . 
We extract Ato^, /+//o the ratio of T{AS) branching 
fractions to B+B~ and B°B°, and the CPT violation 
parameters 5R(cos6') and Q{cos9). Our previous deter- 
mination of Arud using dilepton events from 5.9 fb~^ 
treated /+ //o as a fixed parameter. Otherwise the results 
reported here use the same analysis method and include 
the earlier data and, therefore, supersede the previous 
values. 



Hadronic Event Selection 

Hadronic events are selected from a data set cor- 
responding to 29.4 fb~^ at the T(4S') resonance and 
3.0 fb"^ at an energy 60 MeV below the peak. Hadronic 
events are required to have at least five tracks, an event 
vertex with radial and z-coordinates (where the z axis 
passes through the nominal interaction point, and is an- 
tiparallel to the positron beam) within 1.5 cm and 3.5 cm 
respectively of the nominal interaction point (IP), a total 
reconstructed centre-of-mass (CM) energy greater than 
0.5 W {W is the T(4S') CM energy), a z component of 
the net reconstructed CM momentum less than 0.3 W/c, 
a total ECL energy deposit between 0.025 W and 0.9 W, 
and a ratio i?2 of the second and zeroth Fox- Wolfram 
moments 0] less than 0.7. 

Dilepton Event Selection 



EVENT SELECTION 

The Belle detector, which consists of a silicon vertex 
detector (SVD), a central drift chamber (CDC), aerogel 
Cerenkov counters (ACC), time-of- flight counters (TOP), 
an electromagnetic calorimeter (ECL), and a muon and 
Kl detector (KLM), is described in detail elsewhere 0. 
For electron identification, we use position, cluster en- 
ergy, and shower shape in the ECL, combined with track 
momentum and dE/dx in the CDC and hits in the ACC. 
For muon identification, we extrapolate the CDC track 
to the KLM and compare the measured range and trans- 
verse deviation in the KLM with the expected values. 

Lepton Selection 

The efficiencies for identifying leptons are deter- 
mined from two-photon process data samples: e+e^ 
e^e~e^e~ for electrons and e^e~ — > e^e^/i^/i^ for 
muons. For both cases, the possible degradation of effi- 
ciency due to nearby tracks that are not present in QED 
events but must be considered in hadronic events is ex- 
amined using special hadronic event data samples that 
contain embedded Monte Carlo (MC) lepton tracks. 

We determine the probabilities for misidentifying 
hadrons as leptons using data samples of Ks tt+tt^ 
decays for pions, <j) K^K~ decays for kaons, and 
A — > p7r~ decays for protons. For tracks in the kinematic 
region of the dilepton event selection, the identification 
efficiencies are 92.6% for electrons and 87.0% for muons. 
About 0.1% of pions and kaons, 0.2% of protons and 1.2% 
of antiprotons are misidentified as electrons. About 1% 
of pions and kaons, and 0.2% of protons and antiprotons 
are misidentified as muons. 



Lepton candidates are selected from charged tracks 
that have a distance of closest approach to the run- 
dependent IP less than 0.05 cm radially (drip) and 2.0 cm 
in z {dzyp). At least one r-0 and two z hits are required 
in the SVD. To eliminate electrons from 7 — » e+e~ con- 
versions, electron candidates are paired with all other 
oppositely charged tracks and the invariant mass (assum- 
ing the electron mass hypothesis) Mg+g- is calculated. If 
A/g+g_ < 100 MeV/c^, the candidate electron is rejected. 
If a hadronic event contains more than two lepton can- 
didates, we use the two candidates with the highest CM 
momenta. 

The CM momentum of each lepton is required to be 
in the range 1.1 GeV/c < p* < 2.3 GeV/c. The lower 
cut reduces contributions from secondary (charm) de- 
cay. The upper cut reduces the contribution from non- 
Si? continuum events. The angle of each lepton track 
with respect to the z axis in the laboratory frame must 
satisfy 30° < 6'iab < 135°. This cut rejects tracks 
with large angles of incidence in the SVD and hence 
provides better z vertex resolution. In addition, these 
cuts remove lepton candidates whose particle identifica- 
tion is performed using the endcap KLM or ECL, where 
the performance is not as good as that of the barrel 
sections. Events that contain one or more J/ip can- 
didates are rejected. We calculate the invariant mass 
of each candidate lepton with each oppositely charged 
track (assuming the correct lepton mass hypothesis). 
If the invariant mass falls into the J/tp region, defined 
as -0.15 GeV/c^ < (Mg+g- - Mj/^) < 0.05 GeV/c^, 
-0.05 GeV/c2 < (Af^+^- - Mj/^) < 0.05 GeV/c^, the 
candidate event is rejected. The looser lower cut for 
the electron pair invariant mass is to reject J/ip mesons 
whose calculated mass is low due to bremsstrahlung 
of the daughter electron(s). The opening angle of the 
two leptons in the CM frame is required to satisfy 
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—0.8 < cosO^g < 0.95 to reduce jet-like continuum events, 
SS events where the two reconstructed tracks originate 
from the from the same particle and events with a pri- 
mary lepton and a secondary lepton originating from the 
same B meson. After all cuts are applied, we obtain 
49838 SS and 230881 OS events. The numbers of selected 
dilepton events are summarised in table 

TABLE I: Summary of dilepton events after all cuts. The 
numbers in the off-resonance columns are scaled using the 
luminosity ratio and represent the contribution included in 
the on-resonance data. 



Lepton flavours 


On-rcsonance 


Ofi^-rcsonance 




SS OS 


SS OS 


ee 

MM 
en 


9877 52141 
15503 65435 
24458 113305 


107.4 1513.3 
1464.4 4451.9 
976.3 4403.1 


total 


49838 230881 


2548.1 10368.2 



Az DETERMINATION 

The z-coordinate of each B meson decay vertex is 
determined from the intersection of the lepton track 
with the run-dependent profile of the IP smeared in the 
r-(j) plane by 21 //m to account for the transverse B 
flight length. We define the difference between the z- 
coordinates of the two leptons as Az = z{£~^) — z{l^) for 
OS events and Az = \z{t^) - z{t^) \ for SS events. Az is 
related to the proper-time difference by Az ~ cfi^fAt. 
The Lorentz boost factor of the e"'"e~ CM frame at 
KEKB is /37 = 0.425 0. 

The observed Az distributions have contributions from 
"signal" defined as events where both leptons are pri- 
mary leptons from semileptonic decays of B mesons, and 
"background" where at least one lepton is secondary or 
fake, or the event is from the continuum. The contribu- 
tions from "signal" are theoretically well-defined distri- 
butions convolved with the detector response function, 
which describes the difference between the true Az and 
the measured Az. In order to estimate the detector re- 
sponse function, we employ J/V' decays, where the true 
Az is equal to zero, and whose measured Az distribu- 
tion, after the contribution of background is subtracted, 
yields the response function. 

Candidate J/ijj events are selected using the same cri- 
teria as the dilepton sample except the J/V' veto and the 
cuts on cos 6'!^ are not applied. We define the J/i/) sig- 
nal region as 3.00 GeV/c^ < M^+^- < 3.14 GeV/c^ and 
3.05 GeV/c^ < M^+^- < 3.14 GeV/c^, and the sideband 
region as 3.18 GeV/c^ < M^+^- < 3.50 GeV/c^ for both 
electrons and muons. A linear function is fitted to the 
sideband region of the mass distribution and extrapo- 
lated to the signal region. Using the fit result, we scale 



the Az distribution in the sideband region to the back- 
ground underneath the peak in the signal region, and 
subtract it from the signal region Az distribution. Fig- 
ure Q shows the resulting detector response function. It 
has RMS = 186 /xm in the range |Az| < 1850 ^m. We 
use this histogram as a lookup table in the analysis. 



E 




AZ (iim) 

FIG. 1: Detector response function determined from the J/ip 
data. 



FITTING 

The mixing parameter Am^ and other parameters are 
extracted by simultaneously fitting the Az distributions 
of SS and OS events to the sum of contributions from all 
known signal and background sources. We use a binned 
maximum likelihood method. The background Az dis- 
tributions are obtained from MC and used in the fits in 
the form of lookup tables. 

Signal Az Distributions 

The signal from the neutral B meson pairs originates 
either from B°B° (unm) or from B°B° and B°B° (mix). 
If CPT is assumed to be conserved, the signal distribu- 
tions are expressed by 

|At| 

= N,s fo bl [1 + cos(Am<jAi)] 

= N,s foblefi- '^[l - cosiArud At)]. (4) 
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If the possibility of CPT violation is included, these ex- 
pressions become, 

+1 + Icos^p - 23(cos6') sin(Am<jAi)] 

p-i- = 7V45/o62e--2_^|sin0|2[l-cos(AmdAt)].(5) 

The integrals of P""™ and p™"^ give the time-integrated 
fraction of mixed events 

Isin^px^ 

" \ sm0\^xl + 2 + xl + \cos9\^xj ^ ' 

where Xd = TgoAnid- When CPT is conserved, this be- 
comes the more familiar expression Xd — + 2a;^). 
The values of |cos0|^ and jsin^p are determined from 
5R(cos0) and 3(cos0). 

The signal distribution for charged B meson pairs is 
the same for both the CPT conserving and CPT violat- 
ing cases and given by 

At 

P^hd ^ ^2 chd £_fl 

ZTb+ 

In the equation above, N4S is the total number of T(4S') 
events, fo and /+ are the branching fractions of T{AS) 
to neutral and charged B pairs (/+ + /o — 1), ^0 and 6+ 
are the semileptonic branching fractions for neutral and 
charged B mesons, eu with superscript are the efhciencies 
for selecting dilepton events of charged (chd), unmixed 
(unm), and mixed (mix) origins. The ratio e^^'^: e™™: 
e"^"^ is determined from MC and is fixed in the fit, be- 
cause any detector effect that is not simulated correctly 
should affect events with these origins equally. The Az 
distributions are obtained from these distributions (equa- 
tions ^ and 01 by conversion from At and convolution 
with the emperical resolution function described in the 
previous section. 



the amount of continuum background, we use the off- 
resonance data. This leaves seven parameters to nor- 
malise the fractions of other backgrounds. The first of 
these is Xd which is varied in the fit as given by equa- 
tion 1^ the remaining six are associated with the effi- 
ciencies for selecting these background events which are 
denoted as 

chd chd unm unm mix mix /q\ 
'^SS^ ^OS! ^SS I '^OS ! ^SS I ^os- [°) 

We combine these six types into two: correct-tag (CT) 
which is associated with e'pg, Eq™, £55^, and wrong-tag 
(WT) which is associated with ef^, eg^™, egg'. For both 
CT and WT backgrounds, relative fractions of the three 
background types are fixed according to the MC. We then 
determine the Az distributions for CT and WT back- 
grounds by adding the three corresponding contributions. 

The shapes of the Az background distributions and the 
normalisation of the backgrounds from neutral B events 
depend on Am^. To account for this, we generated two 
samples of generic neutral B MC events, one with Am^ = 
0.469 ps~^ and one with Am^ = 0.522 ps~^. Background 
Az distributions for any value of Am^ are produced by 
linear interpolation between these two MC data sets. 

The MC simulation does not always reproduce the 
hadron showering processes correctly in the kinematic 
region of interest. In order to account for possible dis- 
crepancies, we obtain an overall correction factor by using 
a special control data sample. In this sample events are 
selected in the same way as the dilepton events, except 
that we now require that exactly one lepton passes the 
lepton selection criteria and the other track passes all 
selection criteria except for the lepton identification re- 
quirements. The ratio of data to MC fake rates binned in 
0iab and p is then applied to the control sample to obtain 
the overall fake rates for the dilepton analysis. From this 
method, the hadron misidentification probabilities of the 
MC are increased by 6% for muons and decreased by 5% 
for electrons. 



Background Az Distributions 

The background Az distributions are estimated using 
the MC. A comparison of the Az distribution of J/tp data 
samples between the data and MC shows that the data 
distribution is wider than the MC distribution. A de- 
tailed study showed that after convolving the MC distri- 
bution with a cr = 50 ± 18 /im Gaussian, the distribu- 
tions compared favourably. We smear MC background 
distributions in the same way to compensate for this dis- 
crepancy. 

We categorise the backgrounds into eight types de- 
pending on their sources: charged B pairs, mixed and 
unmixed neutral B pairs, and continuum, each of them 
contributing to both SS and OS events. To normalise 



RESULTS 

In the fit, we fix the parameters r^o = 1.542±0.016 ps, 
and Tg+/rga = 1.083±0.017 0, and we impose the con- 
straint = tb+Itb"- In the fit where CPT is con- 
served, we have a total of five parameters to be fitted. 
They are Am^, f+lfa, two parameters which are related 
to the fractions of CT and WT backgrounds and are ex- 
pressed as TcT = ^'+e«'^/eos and Twt = ^+e«Ves's , 
and an overall normalisation. In the search for CPT vi- 
olation, we add two parameters 5R(cos6') and 3(cos6'). 
The results of the fits are summarised in table Hll 

The values of Ato^ and /+ //q do not vary in a signif- 
icant way when CPT violation is included. We use the 
result from the CPT conserving fit to obtain Anid and 
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TABLE II: Results of fits. 



Fitting method 


CPT conserved 


Allow CPT Violation 


Arud {ps^'-) 


0.503 ± 0.008 


0.503 ± 0.008 


/+//() 


1.01 ± 0.03 


1.02 ± 0.03 


^(cose) 




0.00 ±0.12 


9(cos6') 




0.03 ±0.01 


TcT 


94 ±6 


91 ±5 


TwT 


140 ±2 


139 ±2 




139 {ndf = 86) 


132 {ndf = 84) 



f+ffo- The Az distributions for the SS and OS events 
and the resulting asymmetry are shown in Fig. [21 to- 
gether with the fit results. 
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FIG. 2: Results of the simultaneous fit to the same-sign 
and opposite-sign Az distributions assuming CPT invariance. 
The upper two plots are for the opposite-sign and same-sign 
events, respectively. Signal and background obtained from 
the fit are also shown. The bottom plot is the OS-SS asym- 
metry with the fit result superimposed. 



Systematic errors in the determination of Ato^ and 



Z+Z/o include contributions from uncertainties in the pa- 
rameters that are fixed in the fit: t^+Zt^o, the extra 
50 /zm smearing of the background Az distributions, tbo , 
inclusive D meson branching fractions B — > D'^X and 
B D^X. Also contributing to the systematic errors 
are uncertainties in the continuum contribution deter- 
mined from off-resonance data, lepton misidentification 
probabilities, detector response function, MC statistics, 
background function linear interpolation end points, IP 
profile smearing, the event selection criteria ^lab and dnp, 
and the Az fit range. Uncertainties of the overall z scale 
in the detector and the boost factor (3^ can also con- 
tribute to the systematic error. 

The contribution from the detector response function 
is dominated by the statistics of the J/ip event sample. 
Biases due to approximating the detector response func- 
tion with the Az distribution of J/^ events and the ef- 
fect of the B motion in the T(4S') frame also have some 
contribution. We estimate the uncertainty of the detec- 
tor response function by comparing the results of fits to 
a full MC simulation using three different types of re- 
sponse functions. These are the J/tj} response function 
(extracted in the same way as in data), a true Az re- 
sponse function, and a response function based on the 
proper time difference between the B meson decays (for 
systematic errors associated with the B motion in the 
T(4S') frame). The true Az response function is con- 
structed using signal events and subtracting the true Az 
from the measured Az. The proper time response func- 
tion is constructed in the same manner, except the true 
PjcAt is subtracted. Each of the three types of response 
functions is generated using 35 independent MC data 
sets, resulting in 105 response functions. The fit results 
using each response function are then compared to ex- 
tract the systematic errors. Additionally we vary the 
amount of background subtraction from the J/ip mass 
peak by ±1 a, and repeat the fits. The final detector 
response function systematic error is the combined error 
resulting from the above methods. 

To quantify the systematic error associated with the 
MC statistics, the MC is divided into n = 20 sets and 
the data is fitted using each background distribution Pi, 
{i = l...n). Both the shapes and the efficiency ratios are 
independent for each fit. The fits were repeated with 45 
different pairs of end points for the Anid linear interpo- 
lation. The RMS of the fit results was assigned as the 
error. The systematic error for the IP constraint was es- 
timated by varying the smearing used to represent the 
transverse B fiight length by ±10 /xm. The cuts on the 
variables Az, 6'iab and drjp were varied in the region of 
their default values and the changes in the fit results as- 
signed as errors. The remaining systematic errors were 
calculated by varying the default values by ±1 tr, repeat- 
ing the fits and assigning the differences as errors. The 
contributions from the uncertanies in the z scale and 
were found to be negligible. 
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The errors are summarised in table IIIll 



TABLE III: Systematic errors contributing to the Am^ and 
f+/fo measurements. 

Source 



Amd(ps ^) ,f+//o 



Tb+/tbo (±0.017) ±0.0053 ±0.071 

Detector response function ±0.0047 ±0.021 

Monte Carlo statistics ±0.0036 ±0.011 

50 /im smearing of blig Az (±18 /im) ±0.0032 ±0.015 

Az cut ±0.0030 ±0.010 

Slab cut ±0.0030 ±0.030 

tbo (±0.016 ps) ±0.0012 ±0.002 

driP cut ±0.0010 ±0.003 

Continuum contribution ±0.0008 ±0.001 

Br{B DX){D°: ±4.7%, D+: ±7.9%) ±0.0007 ±0.001 

Fake rate correction (/i: ±3%, e: ±25%) ±0.0007 ±0.002 

Linear interpolation ±0.0005 ±0.020 

IP profile (±10 Aim) ±0.0002 ±0.005 



Quadratic sum 



±0.0097 ±0.086 



To reduce the backgrounds, the analysis was repeated 
with tighter p* cuts, and separate fits were performed 
with ee, efj, and /x^t sub-samples. Deviations from the 
default results were all consistent with statistical fluc- 
tuations. We repeated the fit including the effects of 
Ar/r — 1% and found the shift in results to be negligi- 
ble (w 0.0001 ps-i). 

The largest contributions to the systematic error for 
9(cos0) come from the uncertainty in the extra 50 /xm 
smearing and the ^lab selection criteria and amount to 
±0.03. The systematic error for 5R(cos6') is dominated 
by the uncertainty in the response function. We conser- 
vatively assign an error of ±0.01. 

In summary, we obtain 



Amd 
/+//o 
3?(cos0) 
^(cosS) 



(0.503 ± 0.008(stat) ± O.OlO(sys)) ps^^ 
1.01 ± 0.03(stat) ± 0.09(sys) 
0.00±0.12(stat) ±0.01(sys) 
0.03±0.01(stat) ±0.03(sys). (9) 



Using world avera^s for Anid, the neutral B meson mass 
and decay width 3|, these CPT parameters imply 2] 
the upper limits \mgo — mgo\/mga < 1.16 x 10^^"* and 

ir^o - rsol/rso < o.ii at 90% c.l.. 



SUMMARY 

In summary, we have measured Am^ and the ratio 
of branching fractions for T(4S') decay to B^B~ and 
B^B^ , /+ //o, using inclusive dilepton events. The largest 
contributions to the systematic error for the Am^ mea- 
surement come from the response function and the un- 
certainty in the B meson lifetime ratio T5+ / r^o . Re- 



cent measurements [3 of this lifetime ratio have signifi- 
cantly reduced this systematic contribution since our first 
measurement 0]. The result of this dilepton analysis 
Arud = (0.503 ±0.008 ±0.010) ps'^ is in good agreement 
with the results from other methods used by Belle Q 
and the current world average 0. The error for the 
f+/fo measurement is dominated by the uncertainly in 
tb+ /tb°- We have also obtained new limits on CPT vi- 
olation parameters. 
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